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Using numerical calculations of equi-energy contour plot of the Fourier transform of the spin 
Hamiltonian, we study the magnetic phase diagram (J2 and J3) of the rutile type /3-Mn02, where 
Ji (< 0) is fixed and is the antiferromagnetic interaction along the diagonal direction, J2 is the 
interaction along the c axis, and J3 is the interaction along the a axis. The magnetic phase diagram 
consists of the multricritical point (the intersection J2J3 = J\ and J2 + J3 = 2Ji), the helical order 
along the c axis, the (h = 1/2, k — 0,1 — 1/2) phase, the helical order along the a axis, and the 
phase (h = 0, k = 0, 1 = 1). The shift of the location of the magnetic Bragg peak in the (h, 0, 1) 
reciprocal lattice plane is examined with the change of J2 and J3 in the phase diagram. The shift 
is discontinuous on the first-order phase transition, and is continuous on the second-order phase 
transition. The detail of our magnetic phase diagram is rather different from that reported by 
Yoshimori. 

PACS numbers: 75.10.-b,75.30.Kz 



I. INTRODUCTION 

It is well known that /3-MnC>2 (with the rutile-type 
structure) is one of the systems with a helical spin order 
along the c axis. EricksorpJ found from magnetic neutron 
scattering on /3-Mn02 that the spins are helically ordered 
with a period of 7c along the c axis. The direction of 
spins in the c plane turns from one plane to the next 
plane by an angle of 129° (5tt/7). In 1958, Yoshimori^ 
theoretically demonstrated the origin and stability of the 
helical spin order in /3-Mn02 , based on the experimental 
results from Erickson. The spin Hamiltonian of /3-MnC>2 
is described by a form of the Heisenberg-type, where the 
combination of nearest neighbor (n.n.) exchange inter- 
actions, next nearest neighbor (n.n.n.) exchange inter- 
actions, and so on are included. His success lies in the 
use of the Fourier transform of the spin Hamiltonian in 
the reciprocal lattice space. Then the Fourier transform 
depends only on exchange interactions and wavevectors 
Q. The location of the magnetic Bragg points in the re- 
ciprocal lattice space are uniquely determined from the 
condition that the Fourier transform takes a maximum. 

In spite of the success in theory established by 
Yoshimori P Erickson has not published his detailed data 
of the magnetic neutron scattering of /3-MnC>2 as far as 
we know. Since 1958, there have been several papers on 
the magnetic structure of /3-Mn02 using magneti c n eu- 
tron scatteringP^ an d magnetic x-ray scattering?^ In 
recent years, Sato et al, have reported that the mag- 
netic Bragg peaks appear at the wave vector Q = (1, 0, 
2 + e) in the units of a* (= 2n/a), b* (= 2w/b), and 
c* (= 2ir/c) in the reciprocal lattice space, where e = 
0.297 at 10 K, increases with increasing temperature, and 
reaches 0.2992 just below T N (= 92 K). The value of e 
is rather different from 2/7 derived by Yoshimori P This 
indicates that the helical spin structure along the c axis 
is incommensurate with the c- axis lattice constant. 

Because of the crystal field (distorted octahedron 



formed by 2 ~ ions) in the vicinity of Mn 4+ ion, the 
ground orbital state of Mn 4+ ion (3d 3 , L = 3 and S — 
3/2) is split into the t 2g (de) level (lower energy level, 
triple degenerate) and the e g (d'y) level (upper energy, 
double degenerate). As a result, the ground state is now 
orbital singlet, indicating that the orbital angular mo- 
mentum is quenched. The e g electrons are responsible 
for the metallic conduction, while the localized ti g elec- 
trons are responsible for the magnetism. In recent years, 
Sato et alP have reported the transport properties of a 
single crystal /3-Mn0 2 , such as electrical resistivity, ther- 
mopowcr, Hall effect, and magnetoresistancc. Their data 
show an appreciable anomaly near T/v . This implies that 
there is a strong correlation between conduction e g elec- 
trons and localized ti g magnetic moment through the 
Hund's rule. The DC magnetic susceptibility shows a sig- 
nificant deviation from a molecular-field theory based on 
a localized spin model, which was assumed by YoshimorP 
on his helical spin order. 

In the present paper, we study the magnetic phase dia- 
gram of /3-Mn02 type structure using the molecular field 
theory developed by YoshimoriJ^J where Mn + spins are 
localized. The exchange interactions J\ along the diago- 
nal axis of the system, J2 along the c axis and J3 along the 
a axis are taken into account. Note that J\ is assumed to 
be fixed and be antiferromagnetic. According to Sato et 
al.Pthe intra-atomic t2 g — e g exchange interaction (Hund 
coupling), the transfer interaction between adjacent e g 
orbitals, and the occupancy of the e g orbitals are defined 
by Jiiund, t, and c, respectively. Since | JJ > \cJ Hund \ 
and |J„| > \ct\ for /3-Mn0 2 (n = 1,2, 3)P this means 
that the spin structure of /3-Mn02 below Tjv is not af- 
fected by the effect of cJuund and ct at al. In other words, 
the spin structure is well described by the localized spin 
model with Ji, J2, and J3. We find that our magnetic 
phase diagram ( J2 vs J3) consists of four phases including 
the helical phase along the c axis, the phase with (1/2, 
0, 1/2), the helical phase along the a axis, and the phase 



2 



with (1, 0, 0). The detail of our phase diagram is rather 
different from that proposed by YoshimoriP We use a nu- 
merical calculation approach in finding the distribution 
of the magnetic Bragg peaks in the fixed reciprocal lattice 
planes such as (h, k, I) with one index fixed. To this end 
we calculate the equi-energy contour plot of the negative 
sign of the Fourier transform of the spin Hamiltonian, 
J(h,k,l). The magnetic Bragg peaks are located inside 
the maximum equi-energy contour. The selection rule for 
the location of the magnetic Bragg peaks is the same as 
that derived by YoshimoriP This numerical method has 
an advantage in visualizing the location of the magnetic 
Bragg peaks in the reciprocal lattice space. The nature 
of the phase transitions on the phase boundaries will be 
discussed. 



II. BACKGROUND: GENERAL THEORY FOR 
THE ORDERED SPIN STRUCTURE 

We follow the theory presented by NagamiyaP We con- 
sider a lattice of magnetic atoms such as /3-Mn0 2 . The 
unit cell can be chosen so that it contains one magnetic 
atom. On each magnetic atom, we assume a classical 
spin. Between the spin Si at the position Rj and Sj 
at Rj, there is an Heisenberg-type exchange interaction. 
The Hciscnberg spin Hamiltonian is expressed by 



where 



and 



H — — 2 J (Rjj)Si • Sj, 



R,2j — R^ R j j • 



(1) 



The exchange interaction J(Ry) is not restricted to the 
nearest neighbors. We now use the Fourier transforma- 
tions of the exchange interaction and spin; 



J ( c l) = E ^(Ry)exp(-iq-Ry), 
j= ^S q exp(ig • R 4 ), 



Si = 



(2) 
(3) 



with 



S q = s i exp(-iq • R,), 



where N (~ NiN 2 N 3 ) is the total number of spins, and 
S* = S_ 9 . The position vector Rj is expressed by 



where n\ = 0, l,---,Nx, n 2 — 0,I,---,Y 2 , n 3 = 
0, 1, • • • , Y3, and ai, a 2 , and a3 are the fundamental lat- 
tice vectors. We also define the reciprocal lattice vector 
Gby 

G(h, k, I) = hhi + fcb 2 + lh 3 [h, k, I are integers), 

where t>i, b 2 , and b3 are fundamental reciprocal lattice 
vectors and are given by 

, a 2 x a 3 a 3 x ai a : x a 2 

t>i = 2tt- -,b 2 = 2n- T ,b 3 = 2n-, 



[ai,a 2 ,a 3 ] ' 



[ai,a 2 ,a 3 ] ' 



[ai,a 2 ,a 3 ] 



with 



[ai,a 2 ,a 3 ] = a r (a 2 x a 3 ) = a 2 -(a 3 x a x ) = a 3 -(ai x a 2 ) . 
Noting that 

ai ■ bi = 2-7T, a 2 ■ b 2 = 2ir, a 3 • b 3 = 2ir, 

we have 

G(h, k, I) ■ Hi = 2-K[n\h + n 2 fc + n 3 l) = 2ir x integer. 
The periodic boundary condition for Si leads to 
exp[iq-(iViai)] = I, exp[iq-(Y 2 a 2 )] = I , exp[iq-(Y 3 a 3 )] = f. 
This means that the wavevector q is given by 
q = gibi + q 2 b 2 + g 3 b 3 , 

where 

mi m 2 m 3 

qi = w q2 = w 2 ' q3 = w 2 - 

For convenience we assume that 

Ni Nt N 2 N 2 N 3 N 3 

< mi < — , < m? < — , < m 3 < — , 

2~ ~ 2 2 ~ 2 2 ~ ~ 2 

corresponding to the first Brillouin zone. There are 
N1N2N3 — N wavevectors in the first Brillouin zone. 
The spin Hamiltonian is rewritten as 



i/ = -^J(q)S q -S_ 



(4) 



We look for the lowest minimum of Eq. ([I]) under the con- 
dition that 

s * = s2 = jf E s q ■ s -* ex p ^ - q') • a*]. ( 5 ) 

q,q' 

for any i. Instead of this condition, we impose a milder 
condition 

ns 2 = J2 s ? = 4 E E s q • s - q ' ex pM<i - <*') ■ 



q q' 



R, 



niai + n 2 a 2 + n 3 a 3 (m, n 2 , n 3 are integers) 



— S q • S_ q , 



(0) 
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where we use 

^2 CX P[*(q - q') ' R d = ^^q,q' ■ 

i 

Under this milder condition, the minimum of Eq.Q is 
obtained simply by taking only that q for which J(q) 
has the maximum. Denoting this q by Q (q = -Q being 
equally allowed), we have the minimum value of Eq.Q 
as 



-J(Q)(S Q -S_ Q + S_ Q .S Q ). 



We also obtain 
1 



Si = 



[S Q exp(iQ • Ri) + S_ Q expHQ • R»)]. 



The condition (JsJ) can be written as 

NS 2 = 2S Q • S_ Q + S Q • S Q exp[2iQ • Ri] 
+ S_q • S_q exp[-2«Q • Ri], 

for any R^. This indicates that 

s Q • s Q = 0. 

Here we assume that 

Sq = Rq + Hq, 
S-q = Rq — Hq, 

where Rq and Iq are real vectors. 

Sq ■ S_q = (Rq + Hq) ■ (Rq - Hq) 



(7) 



Rq Rq 



Iq Iq 



Sq • Sq = 



Then we have 



(R q + 2 Iq)-(Rq+zIq) 

Rq • Rq Iq • Iq + 2^i?Q • Iq = 0. (8) 



Rq = ^Q 

Rq Iq = 0. 

Then the vector Rq is perpendicular to the vector Iq, 
and the magnitude Rq is the same as the magnitude Iq . 
From Eq.Q, we get 



R Q 



Iq = \NS 2 . 



The minimum energy E m in is obtained as 

E mm = -NS 2 J(Q). 
The spin vector is expressed by 

;(R Q + ilq) exp(iQ • Ri) 



(9) 



(Rq - ilq) exp(-iQ • Ri)] 
2 [(Rq cos(Q • Ri) - I Q sin(Q • R,)] 



N 

S[(Rq cos(Q • Rj) - Iq sin(Q • R,)] 




FIG. 1: (Color online) The position of Mn 4+ ions (denoted by 
solid circles) and O 2- (denoted by open circles) in /3-Mn02 
(rutile-type structure, space group P4 2 /mmm). a = b = 
4.396A. c = 2.871A Ji, J2, and J3 are the exchange in- 
teractions between Mn 2+ spins. One of the O atom is located 
at (ua,ua), where u = 0.302. 



where Rq and Iq are the unit vectors which are perpen- 
dicular to each other. For convenience, Rq and Iq are 
in the x-y plane. The z axis is perpendicular to the x-y 
plane. Then using the unit vectors e x and e y , we get the 
final result 



Si = 5"[cos(Q • R, + (j))e x + sin(Q • R> 



-y\i 



(10) 



where cj> is the angle between Rq and x axis. 

The spin structure thus derived is the most fundamen- 
tal spin structure and is realized as a result of the mini- 
mum energy state in the classical spin system. The fer- 
romagnetic state (Q = 0) and antiferromagnctic state 
(Q = the zone boundary of the first Brillouin zone) are 
the special case of the spin structures. In general, Q is 
not related to the crystal structure, but is related to the 
details of the exchange interactions. 



III. HELICAL SPIN ORDER FOR /3-Mn0 2 
A. Calculation of J(q) for /3-Mn0 2 

Figure [l] shows the structure (rutile) of /3-Mn02 where 
a = b = 4.3961, c = 2.87ll, and u = 0.302P The ex- 
change interactions J\ , J2 , and J3 are defined by Fig. [l] 
The definition of J 1; J 2 , and J3 are the same as that used 
by Yoshimori. The Neel temperature Tjv is equal to 92 
K. Each Mn 4+ ion and surrounding six 2 ~ ions form a 
cation-occupied deformed octahedron (see Sec. jXJ). 

Here we calculate J(q) for Mn02, where the wavenum- 
ber q, the lattice vectors a(n) (n = 1, 2,. . . , 4) and c are 
given by 



q = {qx,q y ,qz), 

acos[-( 
(0,0,c). 



7T 7T 

a(n) = acos[— (n — 1)], asin[— (n — 1)], 0, 
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Then the expression of J(q) is given by 



Tf . T r c + a(l)+a(2) 1 -c + a(l) + a(2) 1 
J (q) = Ji cxp[zq - ] + Ji cxp[iq ] 

^ 7 r . c + a(2)+a(3) -c + a(2) +a(3) 1 

+ Ji exp[iq • ] + Jx cxp[zq • 



+Ji cxp[iq • 



+Ji exp[iq ■ 



2 




c + a(3) - 


-a(4) 


2 




c + a(4) - 


-a(l) 



] + Ji exp[iq • 
] + Ji exp[iq • 



2 




-c + a(3) - 


ha(4) 


2 




-c + a(4) - 


ha(l) 


2 



+ J 2 exp[iq • c] + J 2 cxp[-zq • c] 

+J 3 cxp[iq • a(l)] + J 3 cxp[-iq • a(l)] + J 3 exp[zq • a(2)] + J 3 exp[-iq • a(2)]. 



(11) 



Then J(q) can be rewritten as 

J(q) = J(h, k, I) — 8 Ji cos(7rft) cos(7rfc) cos(7rZ) 
+ 2 J 2 cos(2ttZ) + 2 J 3 [cos(27rft) + cos(27rfc)] , 

(12) 

where 

= ( — )ft = a h, q y = ( — )fc = a k, q z = { — )/ = c I 
a a c 

with h, k, and I being dimensionless numbers and a* and 
c* being reciprocal lattice constants. 



B. Mathematica programs used in the present 
work 

In order to determine the magnetic Bragg peaks in 
the reciprocal lattice plane, we use the following three 
Mathematica programs, ContourPlot for the equi-energy 
contour plot, FindMaximum for finding maximum, and 
ListVectorPlot3D for drawing the spin directions in the 
lattice points of the real space. 



1. Contour plot program 

This program is used to determine the overview on 
the positions of the magnetic Bragg peaks in the (h, k, I) 
reciprocal lattice plane, where one of h, k, and I are fixed. 
When the values of Ji, J2, and J 3 are given, using the 
Mathematica program [ContourPlot] , we can make a plot 
of the contours of J(h, k,l) = a (constant) in the (h, k, I) 
plane, where the parameter a is changed appropriately 
such that the maximum value of J(h, k, I) is obtained. 
This program is very useful to find the selection rule for 
the position of the magnetic Bragg peaks, as a function 
of Ji, J2, and J 3 . 



2. Finding maximum program 

This program is used to determine the exact position of 
the magnetic Bragg peak in the (h, k, I) reciprocal lattice 
plane, where one of h, k, and I are fixed. For convenience, 
here we assume that I is fixed such that I = 0. Using the 
Mathematica program [FindMaximum] , we find the max- 
imum value of J(h, k, I) for given J 2 and J 3 , where Ji (= 
-1) is fixed, and the regions of h and k are appropriately 
chosen. This program is very convenient when one wants 
to know how the position of the Bragg point changes as 
a function of J2 and J 3 . 

3. ListVectorPlot3D program 

This program is used to draw the spin directions at 
each lattice sites in the real space. The spin structure 
depends on the values of Q. 

IV. SELECTION RULE FOR THE MAGNETIC 
BRAGG PEAKS FROM SIMPLE ANALYSIS 

A. Distribution of magnetic Bragg peaks in the (h, 
k, I = fixed integer) reciprocal lattice plane 

The magnetic Bragg peaks can be located at the 
wavevector q, where J(h, k, I) takes a maximum. When 
I = fixed integer, J(h, k, I) is given by 

J(h,k,l) = 8Ji(-l)' cos(7rft) cos(7r/c) 

+ 2J 2 + 2J 3 [cos(27rft)+cos(27rfc)], 

suggesting that the location of the magnetic Bragg peaks 
in the (h, k, I = fixed integer) reciprocal lattice plane 
depends only on Ji and J 3 ,and is independent of the J 2 . 
The selection rule can be derived as follows, 
(i) The distribution of the magnetic Bragg peaks is the 
same in the reciprocal lattice planes, (h, k, 1), (ft, k, 3), 
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(h, k, 5), (h, k, 7),... 

(ii) The distribution of the magnetic Bragg peaks is the 
same in the reciprocal lattice planes , (h, k, 0), (h, k, 2), 
(h, k, 4), (h, k, 8),... 



B. Magnetic Bragg peaks along the / direction 
with h = fixed integer and k — 

When k = and h is a fixed integer, J(h = fixed 
integer, k = 0, I) is rewritten as 

J(h, 0, 1) = 8(-l) h J! cos(7r/) + 2J 2 cos(27r/) + 2 J 3 , 

which indicates that the location of the magnetic Bragg 
peaks along the (h = fixed integer, k = 0, I) direction 
depends only on J\ and J 2 , and is independent of J 3 . 
The following selection rules can be derived. 

(i) The location of the magnetic Bragg peaks along the 
I direction is the same for (1, 0, I), (3, 0, I), (5, 0, I), (7, 

o,0, •••• 

(ii) The location of the magnetic Bragg peaks along the 
I direction is the same for (0, 0, I), (2, 0, I), (4, 0, I), (8, 

0, 0, •••• 



C. Magnetic Bragg peaks along the / direction 
with h = k — half inetgers 

For h = k = half-integer (= 1/2, 3/2, 5/2, • • •), we 
have 

J(h, h, I) = 2 J 2 cos(2tt/) - 4J 3 , 

leading to the appearance of the magnetic Bragg peaks 
appear at 

I = 1/2, 3/2,5/2,- • • 
for h = k = 1/2,3/2,5/2,7/2,- •• . 



D. Magnetic Bragg peaks along the (0,0,1) 
direction 

For any I, J (0,0, 1) can be expressed by 

J(h = 0, k = 0, 1) = J(9) = 8J1 cos(6>) + 2J 2 cos(26»)+4J 3 . 

with 8 = irl . This indicates that the location of the 
magnetic Bragg peaks along the (0, 0, /) direction, is in- 
dependent of J 3 . The derivative of J(0, 0, 1) with respect 
to 9 is obtained as 



dJ(6) 
d9 



= -4[Ji + J 2 cos(0)] sm(0) = 0. 



Then we get the two solutions. 

(i) Ferromagnetic or antiferromagnetic configurations 
sin (9 = 0, 9 = 0,tt. 



where 

.7(0 = 0) - 8Ji + 2J 2 +4J 3 
J(9 = tt) = -8J1 +2J 2 +4J 3 . 

(ii) Helical spin configuration 

cos^ = — Ji/J 2 ,0 = 9 a = arccos(— Ji/J 2 ), 

under the condition of | Ji/J 2 | < 1, where 

72 

J(9 = 9 ) = -4^-2J 2 +4J 3 . 
The difference is calculated as 

J(9 = 6 ) - J(6 = tt) = -4J 2 (^ - l) 2 , 

•J2 

J(9 = 9 ) - J(9 = 0) = -4J 2 (^ + l) 2 . 

Then the helical spin order appears when 

J(9 = 60) > J(9 = tt), 
J(9 = 9 ) > J(9 = 0). 

If the conditions J 2 < and | Ji/J 2 | < 1 are satisfied, the 
magnetic Bragg peaks appear at the wavevectors denoted 

by 

(h, k, I) - (0, 0, lie), (0, 0, 3±e), (0, 0, 5±e), (0, 0, 7±e), 
where e is defined as 

Ji 



cos(7re) 



E. Magnetic Bragg peaks along the (h, 0, 0) 
direction 



J(h, k = 0, 1 = 0) and its derivative are given by 
J(h, k = 0, 1 = 0) = 8 Ji cos(irh)+2 J 3 cos(2tt/i)+2 J 2 +2 J 3 , 
and 

dJ(h,k = 0,1 = 0) 



dh 



= 8w[Ji — J 3 cos(7r/i)] sin(7r/i), 



respectively. Since |Ji/J 3 | > 1 in Mn0 2 , we have only 
sm(wh) = 0. Since J\ < 0, we have the following selection 
rule. The magnetic Bragg appears only at 

(1,0,0),(3,0,0),(5,0,0),(7,0,0),--- 

along the /i-direction. 
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0.30 
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0.2S 
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FIG. 2: (Color online) Plot of e as a function of J2. J\ = —1 



V. DETERMINATION OF J 2 AND J 3 FROM 
EXPERIMENTAL DATA ON /?-Mn0 2 

A. The ratio J2/J1 

We now calculate the value of e, at which J(h = 1, k = 
0, I = 2 + e) takes a maximum, as a function of J 2 , where 
J 1 = — 1. The value of e is uniquely determined as a 
function J 2 as shown in Fig. [2] When J 2 = —1.6039, we 
find e = 2/7 (Yoshimori 2 ). 

Experimentally, the magnetic Bragg peak is observed 
at (/i = 1, fc = 0, Z = 2.29771) for /3-Mn0 2 . Here we note 
that I = 2.29771 = 2 + e, where e = 2/7 + 0.0120. The 
value of Z is slightly deviated from the value of 2 + 2/7, 
which is predicted by Yoshimori P J(h — 1, k — 0,1) and 
its derivative are given by 

J(h = 1, k = 0, Z) = -8 Ji cos(ttZ) + 2 J 2 cos(2ttZ) + 4J 3 , 

and 

J(/i = l,fc = 0,Z) 



f/Z 



47r[2 Ji sin(7rZ) - J 2 sin(27rZ)] 



respectively. J(/i = 1, k = 0, Z) has a maximum when 



cos(-7rZ) 



J 2 



Z = 0±e,2±e,4±e,6±e, ••• 



where 



cos(7re) 



Ji 
J 2 ' 



Since the magnetic Bragg peak appears at I — 2.29771, 
we have 



Jl 



1.68469, 




FIG. 3: (Color online) Plot of |0| /Tjv as a function of ps, 
where p 2 = 1.68469. 



B. 



Ratio p-j = J3/J1 



The ratio J3/J1 cannot be determined uniquely from 
the experimental results on the magnetic neutron scat- 
tering or x-ray magnetic scattering. The ratio |0| /Tjv 
can be calculated as 



|9| = J(0)_ = 4 + P2 + 2 P3 
Tn J(0o) ^+P2-2p 3 ' 



(13) 



where p 2 = J2/J1, P3 = J3/J1, © is the Curie- Weiss 
temperature (the value is negative for /3-Mn0 2 ) and Tjv is 
the Neel temperature. The value of p 3 can be determined 
from the ratio |0| /Tn since p 2 is already determined. We 
make a plot of |0| /Tn as a function of p$ in Fig.[3| where 
p 2 = 1.68469. 

As far as we know, there have been several reports 
on the experimental values of and Tjv for /?-Mn0 2 ; 
= -316 K, T N = 84 K [Bizette and B. Tsai (1949)], 10 
= -1050 K, T N = 92 K [Ohama and Hamaguchi 
(1971)] E = -783 K, T N = 92 K [Sato et al.^OOO)]^- 
respectively. Using the values of Tjv, and the value 
of P2 (= 1.68469) determined above, the value of p 3 
can be calculated as p$ = 0.537, 0.986, 1.091 using the 
data of Bizette and Tsaif^l Sato et al.p and Ohama and 
Hamaguchip' respectively. 

We note that the boundary of the helical phase is de- 
scribed by (the derivation will be shown later) 



P3 = I/P2 or J 2 J 3 = J\. 



(14) 



indicating that J 2 < since J\ < 0. 



The helical order can exist only when p 3 < 1 /p 2 . When 
p 2 = 1.68469 for Mn0 2 , p 3 should be lower than 0.5935. 
In other words, the ratio |0| /Tjv should be lower than 
4.07853. Since T N = 92 K, this means that |0| should 
be smaller than 375.2 K. The value of |0| by Bizette 
and Tsai seems to be reasonable, while the values of |0| 
obtained by Sato et alP and Ohama and HamaguchP are 
much higher than 375.2 K. Note that Yoshimori^ predicts 
that the boundary is described by p 3 p 2 — 1/2. This 
expression is not correct according to our calculation. 
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FIG. 4: (Color online) (a)(b) Distribution of magnetic Bragg 
points in the (h,k = 0,1) reciprocal lattice plane. Ji = — 1, 
J 2 = -1.68469. J 3 = -0.2. p 3 = J3/J1 = 0.2. (c)(d) 
Distribution of magnetic Bragg points in the (h, k = 0,1) re- 
ciprocal lattice plane. J\ — — 1, J2 = —1.68469. J3 = —0.3. 
P3 = J3/J1 = 0.3. (e)(f) Distribution of magnetic Bragg 
points in the (h, k = 0,1) reciprocal plane around h = 1 
and Z = 2. Ji = -1, J 2 = -1.68469. J 3 = -0.4. 
P3 = J3/J1 = 0.4. 



VI. NUMERICAL CALCULATION 

A. The contour plot of J(h,0,l) in the (ft, 0, Z) 
reciprocal lattice plane 

Using the numerical calculation of the contour plot us- 
ing the Mathematica, we find the location of the magnetic 
Bragg peaks in the (h, 0, 1) reciprocal lattice, where 

J(h, 0, 1) — 8J1 cos(7i7i) cos(7rZ) 

+ 2J 2 cos(27rZ) + 2J 3 [cos(27r/i) + l], 

takes a maximum, where J\ = —1, J 2 = P2J1 = 
— 1.68469. Figures |4|a)-(f) show the contour plot of 
J(h, 0, 1) in the reciprocal lattice plane of (h, 0, 1), where 
■^3 = P3J1 (P3 < 0.5935) and J3 is changed as a param- 
eter (J 3 = —0.2,-0.3,-0.4). As is predicted above, we 
find that the magnetic Bragg peaks appear at 

(0, 0, 1 ± e), (0, 0, 3 ± e), (0, 0, 5 ± e), (0, 0, 7 ± e), . . . , 
(1,0,0 ±e), (1,0,2 ±e),(l,0,4±e), (1,0,6 ±e),..., 
(2,0,l±e),(2,0,3±e),(2,0,5±e),(2,0,7±e),..., 



where e = 0.29771 and is independent of pz {pa, < 
0.5935). Note that the essential results are independent 
of the choice of pz at least for 0.2 < p 3 < 0.4. 



B. The contour plot of J(h,k — h,l) in the 
(h, k = h,l) reciprocal lattice plane 

Using the contour plot of the J(h, k — h, I), we find the 
location of the magnetic Braggs in the (h, h, I) reciprocal 
lattice, where 

J(h, h, I) = 8J1 cos 2 (7r/i) cos(7rZ)+2J 2 cos(27rZ)+4J 3 cos(27i7i). 
The location of the magnetic Bragg peaks is given by 



(0,0,l±e),(0,0,3±e),(0,0,5±e),(0,0,7±e),..., 
(1/2, 1/2, 1/2), (1/2, 1/2, 3/2), (1/2, 1/2, 5/2), (1/2, 1/2, 7/2), . . . , 
(l,l,l±e),(l,l,3±e),(l,l,5±e),(l,l,7±e),..., 
(3/2, 3/2, l/2)(3/2, 3/2, 3/2), (3/2, 3/2, 5/2), (3/2, 3/2, 7/2), . . . , 



in the reciprocal lattice plane of (h, k = h, I), where e = C. The contour plot in the (Zi, k, I) plane with Z = 2n 
0.29771. 

Using the contour plot of the Mathematica, we find 
the location of the magnetic Braggs in the (h, k, I = 2n) 
reciprocal lattice, where 



J(h, k, I = 2n) = 8 Ji cos(7rZi) cos(7rfc) 




-2-10 1 2-2-10 1 



k k 

FIG. 5: (Color online) (a) Magnetic Bragg points in the 
(h, k, I — 2n) reciprocal lattice plane, (b) Magnetic Bragg 
points in the (h, k,l = 2n + 1) reciprocal lattice plane. 
Ji = -1, J 2 = -1.68469. J 3 = -0.3. p 3 = 0.3. 

+ 2J 3 [cos(2tt/i) + cos(2vrfc)] + 2J 2 , 

which does not depend on the index n. Figures [5^a) 
show the location of magnetic Bragg peaks in the in- 
plane (h, k, I = 2n) reciprocal lattice plane. We find that 
the magnetic Bragg peaks appear for h + k = even. This 
selection rule is true for I = 2n with any integer n. 



the maximum of 

J(h = 0,k = Q,l) = 4J 3 + 8 Ji cos(Ztt) 
+ 2J 2 cos(2Z7r), 

by taking the derivative of J(h = 0, k = 0, 1) with respect 
to I. The condition of the local maximum is given by 

cos(7W) = — — -, 
•h 

where |Ji/J 2 | < 1. Then the maximum of J(h = 0, k = 
0,1) is given by 

Jm^(h = 0,k = 0,1) = 4J 3 - 4^ - 2J 2 , (15) 

for the helical order with (h = 0, k — 0, and I) where 
I = 0.70228, 



2. Helical order along the a axis 



D. The contour plot in the (h,k,l) plane with 
I = 2n + 1 

Using the contour plot of the Mathematica, we find the 
location of the magnetic Braggs in the (h, k, I = 2n + 1) 
reciprocal lattice, where 

J(h,k,l = 2n + 1) = — 8 Ji cos(7r/i) cos(7rfc)) 
+ 2J 3 [cos(27r/i) + cos(27rfc)] + 2J 2 , 

which does not depend on the index n. Figure [5jb) show 
the location of the magnetic Bragg peaks in the in-plane 
(h, k, I = In + 1) reciprocal lattice plane, where n is the 
integer. The magnetic Bragg peaks appear for h + k = 
odd for I = 2n+ 1. 

VII. PHASE DIAGRAM IN THE (J 2 , J 3 ) PLANE 

The phase diagram for J\ = — 1 can be determined us- 
ing the programs of the ContourPlot and FindMaximum. 
The phase diagram consists of the four phases; (i) the he- 
lical order along the c axis, (ii) the helical order along the 
a axis, (iii) the ordered phase with h = 1/2, k — 0, and 
I = 1/2, and (iv) the ordered phase with h = 0, k = 0, 
I = 1 is shown in the phase diagram Fig. [6] 

A. Ordered phases 

1, Helical order along the c axis 

The spin vectors in the same ab plane are parallel, i.e., 
h = and k = 0. They screw along the c axis. We find 



The helical spin order where the spin vector in the 
same be plane are parallel and they screw along the a 
axis. We find the maximum of 

J(h, k = 0, 1 = 1) = 2J 2 +2J 3 +8Ji cos(7r/i)+2J 3 cos(27r/i), 

by taking the derivative of J(h, k = 0, 1 = 0) with respect 
to h. The condition of the local maximum is given by 

cos(nh) = — — 

where | Ji/J 3 | < 1. The maximum value is 

Jmax(M = 0,Z = l) = 2J 2 -4^. (16) 

^3 



3. The ordered phase with h = 0, k = 0, and I — 1 

J(h = 0, k = 0, 1 = 1) is given by 

J(h = 0, k = 0, 1 = 1) = -8J1 + 2J 2 + 4J 3 . (17) 



4. The ordered phase with h = 1/2, k = 0, and 1 = 1/2 

This phase corresponds to the MnF 2 -type antiferro- 
magnetic structure (named by Yoshimori 2 ). J(h 
12. k = 0,1 = 1/2) is evaluated as 

J(h = 1/2, k = 0,1 = 1/2) = -2J 2 . (18) 
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B. Boundaries between ordered phases 

1. The phase boundary between the spin order with the 
phase (h — 1/2, k = 0, and 1 — 1/2) and the helical spin 
order along the c axis 

The difference of J(h, k, I) between the helical spin or- 
der along the c axis and the phase (h = 1/2, k = 0, 1 = 
1/2) is given by 

J max (/i = 0,fc = 0, - J{h = 1/2, k = 0, 1 = 1/2) 

J 2 

= 4J 2 -4^. 

So the helical ordered phase is energetically favorable 
when 



J 2 > 



Js' 



with the condition of |</i/J 2 | < 1- When J 3 < 0, this 
inequality can be rewritten as 



J 2 J 3 < J\ 1 
P2P3 < 1. 



(19) 



H - c 



(*=i/2, *=o, i=ia) \\ 



J, =-1 



J-. 



H - a 



FIG. 6: (Color online) Scans A, B, and C in the phase diagram 
(J2, J3) to determine the location of the magnetic Bragg peaks 
in the (h,k,l) reciprocal space. Ji = —1. Note that J2 = 
-1.68469 and J 3 = -0.537 for /3-Mn0 2 (denoted by an open 
circle). 



Then the helical order is stable for J 3 < —1. This in- 
equality is satisfied for the condition of the helical order 

jJi/Jal < 1. 



2. The phase boundary between the spin order with (h = 0, 
k — 0, 1 = 1) and the spin order with (h = 1/2, k — 0, 
I = 1/2; 

The difference of J(h,k,l) between the (h = 0,k — 
0,1 = 1) and the phase (h = 1/2, k = 0, 1 = 1/2) is given 
by 

J{h = 1/2, k = 0, / = 1/2) - J(/i = 0, k = 0, Z = 1) 

= -4J 2 -4J3 + 8J1. 

When 

2J 1 >J 2 + J 3 , (20) 

the phase (ft, = 1/2, k — 0,1 — 1/2) is energetically favor- 
able. When 

2Ji < J 2 + J 3 . (21) 
the phase h = 0, k = 0, Z = l)is energetically favorable. 



C. Multicritical point (J2 = —1 and J3 = —1) and 
phase transition lines 

All four ordered phases merge only at the multicritical 
point at (J 2 = —1 and J 3 = —1), where the line given 
by 2Ji = J 2 + J 3 and the curve J 2 J 3 = J\ intersect with 
other (see the phase diagram of J 3 vs J 2 ). Note that 
J\ = — 1. As will be shown later, the line (2 J\ = J 2 + J 3 ) 
is of the first-order. The lines denoted byJ 2 = — 1 and 
J% = — 1 are of the second-order. 

(1) Using the program of finding maximum in J(q) for 
each point (J 2 , J 3 ) in the phase digram, we find the local 
maximum point (h, 0, 1) in the h-l plane. The location 
of these points is plotted as a function of J 2 when J 3 is 
fixed and as a function of J 3 when J 2 is fixed. 

(2) Using the program of finding maximum in J(q) for 
each point (J 2 , J 3 ) in the phase diagram, we find the 
local maximum point (h,k,0) in the h-k plane. These 
points are plotted as a function of J 2 when J 3 is fixed 
and as a function of J 3 when J 2 is fixed. 



3. The phase boundary between the helical order along the 
a axis and the phase (h — k = 1/2, and 1 = 1) 

The difference of J(h,k,l) between the (h = 0, k = 
0,1 = 1) and the helical order along the a axis is given 
by 

■WO, k = 0, 1 = 1) - J(h = 0, k = 0, 1 = 1) 

= -4J 3 (^-l) 2 . 



1. Scan A 

In the Scan A (see Fig. [6]), we choose J 3 = —0.4 
and J\ = — 1. J 2 is changed as a parameter between 
—4.0 and 0. We make a plot of I vs J 2 in Fig. IWa). 
The discontinuity in I vs J 2 occurs at J 2 = —2.5 [the 
phase boundary (first order) between the phase with 
(h = 1/2, k = 0,1 = 1/2) and the helical order along the c 
axis], where the relation J 2 J 3 = Jf is satisfied (Ji = — 1). 
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FIG. 8: (Color online) h vs J 3 along the scan B (—4 < J 3 < 
and J2 = 1.5). Ji = — 1. k = and I = 1. 



(h = 0, fc = 0, I = 1)]. We find that the relation of h vs 
J3 is well described by 

Ji 



cos(wh) = — 



■h 



(23) 



FIG. 7: (Color online) (a) / vs J2 and (b) ft vs J2 along the 
scan A (-4 < J 2 < and J 3 = -0.4). fc = 0. Ji = -1. 



The value of /i is independent of J2. For example, when 
h = 0.2, we have J3 = —1.236. The plot of I vs J2 (which 
is not shown here) indicates that I remains constant (I = 
1) for -4 < J 3 < 0, where J 2 = 1.5. 



3. Scan C 



We find that the relation of I vs J2 is well described by 



cos(7r£) 



Jl 

J 2 ' 



(22) 



for —2.5 < J2 < — 1 (the helical phase along the c axis). 
For example, when I = 0.8, we have J2 = —1.24 as shown 
in Fig. [7]ja). We note that I is equal to 1 at J2 = —1 [the 
phase boundary(second order) between the helical phase 
along the c axis and the phase with [h — 0, k = 0, 1 = 1)] 
and remain unchanged (I = 1) for —1 < J2 < 0. Wc 
also make a plot of h vs J2 in Fig. ^b) ■ The value of h 
undergoes a sudden change from h = 0.5 to at J2 = 
—2.5 (at the phase boundary between the phase with 
(h = 1/2, k = 0, 1 = 1/2) and the helical order along the 
c axis) and remains constant (h = 0) for —2.5 < J2 < 0. 
Note that there is no change of h at J 2 = — 1- 



In the Scan C, we choose J 2 = —0.4 and J\ = — 1. 
J3 is changed as a parameter between -3.5 and 0. We 
make a plot of h vs J3 and I vs J3 in Figs. |9ja) and 
(b). As shown in Fig. |9^a), the value of h decreases with 
increasing J3, showing an abrupt decrease at J3 = —2.5 
(which is denoted by J2J3 = Jf ) [the phase boundary 
(first-order) between the phase with (h = 1/2, fc = 0, / = 
1/2) and the helical order along the a axis]. The value of 



h decreases with further increasing J 3 , following Eq. ( 23 1 , 
and reduces to zero at J 3 = —1 [the phase boundary 
(second-order) between the helical order along the a axis 
and the phase with [h = 0, fc = 0, 1 = 1)]. 

In Fig. gb), the value of I increases with increasing 
J 3 . It shows an abrupt change from i=l/2toZ = lat 
J3 = —2.5 at the phase boundary (first-order) between 
the phase with (h = 1/2, fc = 0, 1 = 1/2) and the helical 
state along the a axis. 



2. Scan B 

In the Scan B (see Fig. [6]), we choose J 2 = 1.5 and 
Ji = — 1. J3 is changed as a parameter between -4.0 
and 0. We make a plot of h vs J2 in Fig. |8j The value 
of h decreases with increasing J3 and reduces to zero at 
J3 = — 1 [the phase boundary (second-order) between 
the helical order along the a axis and the phase with 



VIII. CONTOUR PLOT OF (h, 0,1) FOR 
TYPICAL POINTS IN THE (J 2 , J 3 ) PHASE 
DIAGRAM 

Here we show the contour plot of J(h, 0, /) for the typ- 
ical points in the (J2, J3) phase diagram (see Fig. [To| . 
The maximum points in the countour plot of J(h, 0, 1)) 
correspond to the magnetic Bragg peaks in the (h, 0, 1) 
plane. 
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FIG. 9: (Color online) (a) h vs J3 and (b) I vs J3 along the 
scan C (-3.0 < J 3 < and J 2 = -0.4). Ji = -1. k = 0. 
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FIG. 11: (Color online) Distribution of Bragg peaks in the 
(h,0,l) reciprocal lattice plane. Ji = —1. J2 = —0.5, and 
J3 is changed as a parameter, (a) Point A (J3 = —0.5), (b) 
Point B (J 3 = -1.5), (c) Point C (J 3 = -2.0), (d) Point D 
(J 3 = -3.0), and (e) Point E (J 3 = -4.0). 



A. The points A, B, C, D, and E with Ji = 1.0 

and J2 = —0.5 

Typical contour plot of (h,0,l) at the points A, B, C, 
D, and E in the (J2, J3) phase diagram are shown in 
Fig. 11 The point A ( J 3 = —0.5) is in the phase with (h 
= 0, k = 0, 1 = 1). The point B ( J 3 = -1.5) is in the helical 
phase with the a axis. The point C ( J3 = —2) is the phase 



boundary between the phase with (h = 0, k = 0, / = 1) 
and the helical phase with the a axis. The points D 
(J3 = —3.0) and E (J3 = —4.0) are in the phase with 
(h=l/2,k = 0,l = l/2). 



FIG. 10: (Color online) Phase diagram of J3 vs J2: Ji = 
— l.H-a: helical order along the a axis. H-c: helical order along 
the c axis. The ordered phase with (h — 1/2, k = 0, 1 = 
1/2). The ordered phase with (ft = 0, k = 0, / = 1). The 
multicritical point is at (J2 = —1, J3 = —1). The calculations 
of the contour plots are made at the points (A, B, C, • • •, P). 
The line J2J3 = J? is the first-order phase boundary. The 
line (J2 = —1, J 3 > —1) and the line (J3 = —1, J2 > —1) are 
the second-order phase boundary. 



B. The points F and G with Ji = -1.0 and J 2 = 1.5. 

Typical contour plot of (h, 0, 1) at the points F and G 
in the (J2, J3) phase diagram are shown in Fig. 12 The 
point F ( J3 = —0.5) is in the phase with [h = 0, k = 
0, 1 — 1). The point G (J 3 = —1.5) is in the helical phase 
with the a axis. 
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FIG. 12: (Color online) Distribution of Bragg peaks in the 
(ft, 0, Z) reciprocal lattice plane. Ji = —1. J2 = 1.5, and J3 
is changed as a parameter, (a) F (J3 = —0.5) and (b) G 
(J 3 = -1.5). 
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FIG. 13: (Color online) Distribution of Bragg peaks in the 
(h,0,l) reciprocal lattice plane. Ji = —1. J2 = —1.5, and J3 
is changed as a parameter, (a) H (J3 = —0.2), (b) I (J3 = 
-0.5), (c) J (J 3 = -0.667), (d) K (J 3 = -1.2), and (e) L 
(J 3 = -1.5). 



C. The points H, I, J, K, and L with Ji = —1.0 and 
J2 = -1.5 

Typical contour plot of (ft, 0, 1) at the points H, I, J, 
K, and K in the (J2, J3) phase diagram are shown in 
The point H (J 3 = —0.2) is in the helical phase 
ihe c axis. The point I ( J3 = —0.5) is in the helical 
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FIG. 14: (Color online) Distribution of Bragg peaks in the 
(ft, 0, 1) reciprocal lattice plane. Ji = —1. J2 = —2, and 
J3 is changed as a parameter, (a) M (J3 = —0.2), (b) N 
(J 3 = -0.5), (c) O (J 3 = -0.7), and (d) P (J 3 = -0.8). 



phase with the c axis and is on the line (J2 + J3 = 2Ji). 
The point J (J 3 = —0.667) is on the phase boundary 
between the helical phase with the c axis and the phase 
with (ft = 1.2, k = 0,1 = 1/2). The points K (J 3 = -1.2) 
and L (J3 = —1.5) are in the phase with (ft = 1/2, k = 
0,1 = 1/2). 

D. The points M, N, O, and P with Ji = -1.0 and 

J 2 = -2 

Typical contour plot of (ft, 0, /) at the points M, N, O, 
and P in the ( J2, J3) phase diagram are shown in Fig. 14 
The point M (J3 = —0.2) is in the helical phase along 
the c axis. The point N (J3 — —0.5) is on the phase 
boundary ( J2 + J3 — 2Ji) between the helical phase with 
the c axis and the phase with (ft = 1.2, k = 0,1 = 1/2). 
The points O (J 3 = -0.7) and P ( J 3 = -0.8) are in the 
phase with (ft = 1/2, k = 0,1 = 1/2). 



IX. 3D SPIN STRUCTURES 

What is the three dimensional (3D) spin structure 
which is characterized with the wavevector of the mag- 
netic Bragg peaks? The vector of spin at the site of 
the real lattice space is given by 

Si = 5[cos(Q • Ri)e x + sin(Q • Hi(f>)e y ] 

where S = 3/2, we assume that the phase factor cf> is 
equal to zero, and Q is defined as 

Q = (ha*,ka*,lc*). 
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FIG. 15: (Color online) structure in the phase (h — 0, fc 
0,1 = 1). 




FIG. 16: (Color online) Spin structure in the phase (h 
1/2, fc = 0,; = 1/2). 



A. Phase with (h = 0, fc = 0, Z = 1) 



FIG. 17: (Color online) Incommensurate spin structure for 
the helical order along the c axis with the phase (0, 0, 0.7023). 




FIG. 18: (Color online) Commensurate spin structure for the 
helical order along the c axis, with the phase (0, 0, 1 = 5/7), 
which is proposed by YoshimoriP 



A typical spin structure for the phase with (0, 0, 1) is 
described in Fig. [15] 



D. Helical ordered phase along the c axis 
(Yoshimori); commensurate structure 



B. Phase with (h = 1/2, fc = 0,1 = 1/2) 

A typical spin structure for the phase with (1/2, 0, 
1/2) is described in Fig. 16 



A typical spin structure with I = 0.71421 is commensu- 
rate with the lattice structure and described by Fig. [18] 
The magnetic unit cell along the c axis is seven times 
larger than the unit cell of the crystal unit cell. 



C. Helical ordered phase along the c axis; (0,0,1) 



E. Helical ordered phase along the a axis: (h, 0, 1) 



The value I is obtained as I = 0.7023 for J 2 = -1.68469 
from Eq. ( 22 ) . A typical spin structure is incommensu- 
rate with the lattice structure and described by Fig. 17 



When J3 = —2.0 and J\ = — 1, ft. is equal to 2/3 from 
Eq. 23 1. Then a typical spin structure is described by 

Fig. m 
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FIG. 19: (Color online) Spin structure for the helical order 
along the a axis with (h ma x = 2/3, 0, 1). J\ = — 1, and 
,h = -2.0. 



X. DISCUSSION: NATURE OF EXCHANGE 
INTERACTIONS IN THE RUTILE-TYPE /3-Mn0 2 

A. Overview 

We show that the magnetic phase diagram with Ji = 
— 1 is uniquely determined by the combination of ,J 2 and 
J3, where the sign of J3 is mainly negative and the sign 
of J 2 is changed between negative to positive. For fj- 
MnC>2, we have p 2 = 1.68469 and as P3 = 0.537, leading 
to J 2 = -1.68469 and J 3 = -0.537 when J x = -1. This 
point ( J2, J3) for /3-Mn02 is in the helical order along the 
c axis, as shown in the magnetic phase diagram (Fig. [6]). 
In Fig.[20|[a), we show the structure of /?-Mn02- Figure 



20 



b) shows the distorted octahedron where one Mn 4+ 
ion (cation) at the point Bi, is surrounded by six O 2- 
ions (anion) at the points Oi, O2, O3, O4, O5, and 06- 
The cation-cation separation (A1C1, BiD^ along a [001] 
axis is considerably smaller (2. 871 A). Because of the 
distorted octahedron formed by O 2- ions in the vicin- 
ity of Mn 4+ ion, the ground orbital state of Mn 4+ ion 
(3d 3 , L — 3 and S — 3/2) is split into the t 2g (de) level 
[lower energy level, triple degenerate; d(xy), d(yz), and 
d(zx) states] and the e g (dj) level [upper energy, double 



degenerate; d(3z 2 



2 ), d(x 2 - y 2 ))]. The t 2g state is 



occupied by three electrons with spin up state | +). The 
t g state is empty. As a result, the ground state is orbital 
singlet, indicating that the orbital angular momentum is 
quenched. In /3-Mn02, the e g electrons are responsible 
for the metallic conduction, while the localized t 2q elec- 
trons are responsible for the magnetism. In Fig. |20^b) , 



the (3z 2 



orbital axis at the point Bi coincides with 




5 Bi0 6 (the z axis). The d(x 2 - y z ) axes coincide with 
C^BiC^ (the x axis) and the OiBiC^ (the y axis). The x 
and y axes are perpendicular to the z axis. However, the 



FIG. 20: (Color online) (a) Crystal structure of /3-Mn0 2 . 
ZA1O4B1 = ZC1O1B1 = 130.62°. ZA1O4A2 = 111.60°. 
ZB1O3D1 = 98.77°. The distance (in the units of A) 
is denoted by numerical value. (b) Distorted octahedra 
of ^-Mn0 2 . ZO1B1O4 = 98.77°, ZO1B1O3 = 81.225°, 
Z0 5 0i0 6 = 89.590°, ZOi0 6 5 = 45.205°, and ZO1B1O5 = 
ZO4B1O5 = 90.0°. The distance (in the units of A) is denoted 
by numerical value. 



x axis is not perpendicular to the y axis (98.78°, 81.22° 



B. Origin of the direct exchange interaction J2 

We consider the exchange interaction J 2 between the 
points Bi and Di (or between the points Ai an d C i ). 
The distance BxDj (Aid) is 2.87lA. In Fig. |(Jb), 
one d(xy) orbital (<2g) from the point Bi bisects the 
angle ZO1B1O3 and meets at the middle point of the 
edge O1O3. The other d(xy) orbital from the point Di 
(the center of neighboring octahedron) bisects the angle 
ZO1D1O3 and meets at the middle point of the edge 
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O1O3. According to GoodenoughJ^ the predominant 
interactions between neighboring cations whose cation- 
ccupied octahedra share an edge, are assumed to be di- 
rect cation-cation exchange interaction. Thus the inter- 
action J 2 is antiferromagnetic (Heitler-London type) for 
/3-Mn0 2 . 



C. Origin of superexchange interaction J\ 



We consider the interaction between Mn 4+ at the point 



Bi and the Mn 4+ at the point Ci as shown in Figs. [20 
(a) and (b), where the distance C1B1 is 3.424A Accord- 
ing to Goodenoughp-Q this interaction (Ji) is a superex- 
change one since the cation-occupied octahedra share a 
common corner (point Oi). The angle C1O1B1 is equal 
to 130.62°, this interaction is antiferromagnetic for /3- 
Mn0 2 . According to Goodenoughf^ when the cation- 
occupied octahedral share a common corner, there can 
be no direct overlap of neighboring cation orbitals and 
therefore there is no cation-cation interactions. 

What is the origin of Ji? The point Ci is the center 
of the neighboring octahedron. The direction of C1O1 
is the z axis of this octahedron. A p orbital of Oi is 
expected to be directed toward Ci (p a orbital) so as to 
overlap the d(3z 2 — r 2 ) orbital of the point Ci, where the 
p orbital (p x , p y , p z states) is called as p a orbital when 
the principal axis of the p orbital coincides with the di- 
rection of the bond. A partial covalent bond between 
the d(3z 2 — r 2 ) orbital and the p a orbital can be formed. 
Then the charge transfer occurs from the p a orbital with 
the spin- up state |f) to the e g state. Consequently, the 
spin of Mn 4+ at the point Ci is still in the spin-up state, 
while the resulting spin of 2 ~ at the point Oi is in the 
spin-down state. The remaining p a state at the point Oi 
is magnetically coupled with the ti g state of the Mn 4+ 
at the point Bi. When this coupling (denoted as Kq) 
is antiferromagnetic, then the suprexchange interaction 
between Mn 4+ at the point Ci and Mn 4+ at the point Bi 
is ferromagnetic. On the other hand, when this coupling 
Kq is ferromagnetic, then the superexchange interaction 
between Mn 4+ at the point Ci and Mn 4+ at the point 
Bi is antiferromagnetic. Here we note that the angle 
a = ZC1O1B1 is equal to 130.62° for /3-Mn0 2 , which is 
very different from 90°. If a = 90°, the p a orbital on 
the bond O1C1 coincides with the p^ orbital on the bond 
CqBi (the y axis), where the p orbital {p x , p y , p z states) 
is called as p„ orbital when the principal axis of the p or- 
bital is perpendicular to the direction of the bond. The 
Ptt orbital is coupled with the i 2s state of the Mn 4+ at 
the point Bi, since the e g state is empty. According to 
the Goodenough-Kamamori- Anderson rulep^^ the in- 
teraction Kq is antiferromagnetic, since the p^ orbital is 
not orthogonal to the t2 g orbital. If a = 180°, the p a 
orbital is coupled with the t2 g orbital of the Mn 4+ at the 
point Bi the interaction Kq is ferromagnetic, since the 



Pu orbital is orthogonal to the t2 g orbital. 

The sign of Kq is dependent on the value of a. There 
may be a critical angle a c . Kq is ferromagnetic for a > a c 
and Kq is antiferromagnetic for a < a c . Experimentally, 
the superexchange interaction J\ is antiferromagnetic for 
/3-Mn0 2 , which means that Kq is ferromagnetic. The 
critical angle a c is lower than 130.62°. 

The direct cation-cation exchange interaction (J 2 ) is 
expected to be stronger than the superexchange in- 
teractions J\. The competition between J\ and J 2 
can lead to a complicated compromise magnetic order; 
J2/J1 = 1.68469 for /3-Mn0 2 . Note that the discussion 
of OsmoncP^ on the nature of J2 may be inappropriate. 



D. Origin of superexchange interaction J3 

We consider the interaction between Mn 4+ at the point 
Ai and the Mn 4+ at the point A 2 , as shown in Figs. 20 'a) 
and (b), where the dista nce AiA 2 is 4.396A. Accord- 
ing to Goodenough j 11 1 12 1 this interaction (J 3 ) is a su- 
perexchange one (cation-anion-cation) since the cation- 
occupied octahedra share a common corner (point O4). 
The angle ZAi04A 2 is equal to a = 111.60°, the dis- 
tance 04A 2 is 3.343A, and the distance O4A1 is 1.878A 
This interaction is experimentally antiferromagnetic for 
/3-Mn0 2 ; p% = J3/J1 = 0.537. This means that the crit- 
ical angle a c is between 90° and 111.60°. 



XI. CONCLUSION 

We have studied the phase diagram of ( J 2 vs J3) with 
J 1 = —1 in the rutile type /3-Mn0 2 by using the equi- 
energy contour plot. The distribution of the magnetic 
Bragg peaks can be clearly visualized. The magnetic 
phase diagram consists of the multricritical point (the 
intersection J2J3 — J 2 and J 2 + J 3 = 2J X ), the helical 
order along the c axis, the (h = 1/2, k = 0,1 = 1/2) 
phase, the helical order along the a axis, and the phase 
(h = 0, k = 0, / = 1). The phase transition is of the 
first order between the (h — 1/2, k = 0,1 = 1/2) phase 
and the helical order along the c axis, and between the 
(h = 1/2, k = 0, 1 = 1/2) phase and the helical order 
along the a axis. The phase transition is of the second 
order between the phase (h = 0, k = 0, 1 = 1) and the 
helical order along the c axis, and between the phase 
(h = 0, k = 0, / = 1) and the helical order along the a 
axis. 



Acknowledgments 

We are grateful to Prof. H. Sato for useful discussions 
on the itinerant nature of e g electrons in /3-Mn0 2 . 



16 



* suzu ki@binghamto n.edu| 
t itsuko@binghamton.edu 

1 R.A. Erickson, Phys. Rev. 90, 779 (1953). 

2 A. Yoshimori, J. Phys. Soc. Jpn. 14, 807 (1959). 

3 N. Ohama and Y. Hamaguchi, J. Phys. Soc. Jpn. 30, 1311 
(1971). 

4 M. Regulski, R.Przenioslo, I. Sosnowska, and J.-U. Hoff- 
mann, Phys. Rev. B 68, 172401 (2003). 

5 M. Regulski, R.Przenioslo, I. Sosnowska, and J.-U. Hoff- 
mann, J. Phys. Soc. Jpn.73, 3444 (2004). 

6 H. Sato, K. Wakiya, T. Enoki, T. Kiyama, Y. Wak- 
abayashi, H. Nakano, and Y. Murakami, J. Phys. Soc. Jpn. 
70, 37 (2001). 

7 H. Sato, Y. Kawamura, T. Ogawa, Y. Murakami, H. 
Ohsumi, M. Mizumaki, and N. Ikeda, Physica B 329 - 
333, 757 (2003). 

8 H. Sato, T. Enoki, M. Isobe, and Y. Ueda, Phys. Rev. B 



61, 3563 (2000). 
9 T. Nagamiya, Solid State Physics, edited by F. Seiz, D. 
Turnbull, and H. Ehrenreich, 20, 305 (Academic Press, 
New York, 1967). 

H. Bizette and B. Tsai, Colloque sur la polarization de 
la matere cited by Lidiard (Paris, C.N.R.S., 1949) 164, 
Reports on Progr. in Phys. 17, 201 (1954). 

1 J.B. Goodenough, Phys. Rev.117, 1142 (1960). 

2 J.B. Goodenough, Magnetism and the Chemical Bond 
(John- Wiley & Sons, New York, 1963). 

3 J. Kanamori, J. Phys. Chem. Solid 10, 87 (1959). 

4 P.W. Anderson, Solid State Physics, edited by F. Seiz 
and D. Turnbull (Academic Press, New York and London, 
1963) 14, 99. 

5 W.P. Osmond, Proc. Phys. Soc. 87, 335 (1966). 



